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Abstract 

We extend the formahsm relating electromagnetic form factors to transverse 
quark charge densities in the light-front frame to the case of a spin-3/2 baryon 
and calculate these transverse densities for the A(1232) isobar using lattice 
QCD. The transverse charge densities for a transversely polarized spin-3/2 par- 
ticle are characterized by monopole, dipole. quadrupole, and octupole patterns 
representing the structure beyond that of a pure point-like spin-3/2 particle. We 
present lattice QCD results for the A-isobar electromagnetic form factors for 
pion masses down to approximatively 350 MeV for three cases: quenched QCD, 
two-degenerate flavors of dynamical Wilson quarks, and three flavors of quarks 
using a mixed action that combines domain wall valence quarks and dynamical 
staggered sea quarks. We extract transverse quark charge densities from these 
lattice results and find that the A is prolately deformed, as indicated by the 
fact that the quadrupole moment G'_e2(0) is larger than the value —3 character- 
izing a point particle and the fact that the transverse charge density in a A^ of 
maximal transverse spin projection is elongated along the axis of the spin. 

Key words: Electromagnetic form factors. Electric and Magnetic Moments, 
A-resonance 

PACS: 13.40.Gp, 13.40.Em, 14.20.Gk 



1. Introduction 

The question of how the structure of mesons and baryons arises from the 
interaction among their quark and gluon constituents is at the forefront of con- 
temporary research in hadron physics. Key observables in this field are the 
electromagnetic (e.m.) form factors (EEs), which yield the distribution of the 
quark charges in a hadron. The nucleon e.m. EEs have been especially thor- 
oughly studied recently at electron-beam facilities, such as Jefferson Lab, MIT- 
Bates, and MAMI, see Refs. [H, [2, Q for latest reviews. 
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In systems like nuclei or atoms, in which the spatial size R is large com- 
pared to the system's Compton wavelength l/M, electromagnetic form factors 
are essentially three dimensional Fourier transforms (or distorted wave analogs 
thereof) of ground state charge density distributions, and thereby provide valu- 
able physical insight into the structure of the ground state of the system in the 
lab frame. In hadrons, however, where R ~ l/M, the Fourier transform argu- 
ment is inapplicable, and there is no known way to relate form factors to charge 
densities in the lab frame. However, when the hadron is viewed from a light 
front, there is a simple and consistent field-theoretic density interpretation of 
FFs as the Fourier transform of the spatial distribution of the quark charge in 
the plane transverse to the line-of-sight. Equivalently, in the infinite momentum 
frame, the quantity that plays the role of M in the lab frame Fourier transform 
argument is the light-front momentum , which goes to infinity, so the form 
factor is just the two-dimensional Fourier transform of the transverse density 
in the infinite momentum frame In this way, the transverse quark charge 
densities have been mapped out in the nucleon 0, Q , deuteron Q , and pion 0] 
based on empirical FFs. It is important to note, however, that physics in the 
lab and infinite momentum frames is significantly different and hence there is no 
simple relation between the two-dimensional transverse densities in the infinite 
momentum (or light front) frame and the three-dimensional density in the lab 
frame, and we will return to this issue in the discussion of deformation. 

In this work we address the e.m. FFs and transverse charge densities of the 
A(1232) isobar, the lightest nucleon excitation. The nucleon-to-A transition 
has been well measured experimentally over a large range of photon virtualities, 
see Ref. 0] for a recent review. It is dominated by a magnetic dipole transition, 
while the electric and Coulomb quadrupole transitions were found to be small 
(in the few percent range as compared to the magnetic dipole transition). These 
measurements have made it possible to quantify the deformation of the N A 
transition charge distribution 

On the other hand, the information on the e.m. FFs of the A itself is scarce. 
Because of the tiny lifetime of the A, it is of course very difficult if not impossible 
to access these quantities directly in experiment. Fortunately, however, as in the 
case of the N ^ A FFs, it is now feasible to calculate the A FFs using lattice 
QCD. First lattice QCD results in the quenched approximation were presented 
in Ref. [lo| . and first results on the A FFs using dynamical quarks were briefly 
reported in [llj . In this more extended publication we provide further details on 
the lattice evaluation of these FFs as well as the interpretation of the calculated 
FFs in terms of the transverse charge densities. 

The outline of the paper is as follows: Section[2l contains a general discussion 
of the e.m. interaction of a spin-3/2 system along with the definition of its 
e.m. moments and FFs. In Section [31 we find the specific ('natural') values for 
electromagnetic moments that an elementary (pointlike) spin-3/2 particle would 
possess. In Section HI we introduce the light-front helicity amplitudes for the 
e.m. vertex of a spin-3/2 system, and express them in terms of the corresponding 
e.m. FFs. In Section[5l we calculate quark transverse charge densities of a spin- 
3/2 system in terms of the light-front helicity amplitudes. We also calculate 
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the values of the electric dipole, quadrupole and octupole moments of these 
transverse charge densities, and show that for a pointlike spin-3/2 particle they 
vanish. In Section [51 wc describe our three different types of lattice calculations 
of the A(1232) cm. FFs and compare their results. In the first type, we use 
Wilson fermions in the quenched approximation. The second type uses two- 
degenerate flavors of dynamical Wilson quarks {Np — 2). The third type uses a 
mixed action combining domain wall valence quarks with dynamical staggered 
sea quarks including light degenerate up and down quarks and heavier strange 
quarks (Np = 2+1). In Section [71 we use these lattice QCD results to extract 
the transverse quark charge densities in the A resonance. We summarize our 
results, discuss their relation to models and present our conclusions in Section[51 
The lattice results for the A e.m. FFs are tabulated in an Appendix. 



2. The 7*AA vertex and form factors 



Let us consider the coupling of a photon to a A, shown in Fig.[T] The matrix 
element of the electromagnetic current operator J'' between spin-3/2 states can 
be decomposed into four multipole transitions: a Coulomb monopole (EO), a 
magnetic dipole (Ml), a Coulomb quadrupole (E2) and a magnetic octupole 
(M3). We first write a Lorcntz-covariant decomposition for the on-shell 7* A A 
vertex, which exhibits manifest electromagnetic gauge-invariance [^, : 



(Ab',A')|J''(0)|A(p,A)) 



a') 



2, g"/ 



2\ rr 



(2Ma)2_ 

ui3{p, A), 



2Ma 



(1) 



where Ma = 1.232 GeV is the A mass, Ua is the Rarita-Schwingcr spinor for 
a spin-3/2 state, and A (A') are the initial (final) A helicities. Furthermore, 
-^1 2 3 4 the 7* A A form factors, and ^i*(0) = ca is the A electric charge in 
units of e (e.g., = +1). For future reference, we also define the quantity 
r^QV(4Mi). 

A physical interpretation of the four electromagnetic A — > A transitions can 
be obtained by performing a multipole decomposition [isl . [T^ . For this purpose 
it is convenient to consider the Breit frame, where p = —p' = —q/2. Further- 
more, we choose q along the z-axis and denote the initial (final) A spin projec- 
tions along the z-axis by s (s'). In this frame, the matrix elements of the charge 
operator define the Coulomb monopole (charge) and Coulomb quadrupole form 
factors as : 



(|,.'|J"(0)| -|,.) 



{2Ma)5,,, {(<5,±3 Gbo(Q') 

,±^-6,±i) G£2(Q2)1.(2) 
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Figure 1: The 7* AA vertex. The four-momenta of the initial (final) A and of the photon are 
given by p (p') and q respectively. The four-vector indices of the initial (final) spin-3/2 fields 
are given by f} {a), and ^ is the four- vector index of the photon field. 



Using Eq. ([T]), wc can express the Coulomb monopole and quadrupole form 
factors in terms of F* 2348-8: 

Geo = {F* -tF;) + 1^tGe2, (3) 

Ge2 = {f:-tf;)-]^{i + t){f;-tf:). (4) 

Analogously, the matrix elements of the current operator define the magnetic 
dipolc and magnetic octupole form factors. For the transverse spherical com- 
ponent J+i = — --i=(J^ -|- ij'^) we obtain : 

(|,.'|J+i(0)| -|,.)^(-y2)(2AfA)^ 

X \J^s'+^^s+^ +K'-i5s-l + Gmi{Q^) 
-^r (<5,,+3 5, + . +<5,,_i <5,_3 - V3(5,, + i <5,_i) GA/3(g')| (5) 

Using Eq. ([1]) , we can express the magnetic dipole and octupole form factors in 
terms of 2 3 4 as : 

Gmi = {F*+f;) + ^tGm^. (6) 

Gm3 = {f* + f;)~]^{i + t){f; + fi). (7) 

At = 0, the multipole form factors define the charge (ba), the magnetic 
dipole moment {(J-a), the electric quadrupole moment (Qa), and the magnetic 
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octupole moment (Oa) as : 



^J■A 
Qa 
Oa 



Geo{0) = F*{0), 

-GMi(0) = -^[eA + F;(0)] , 



2Ma 

e 
M 



2 Gb2(0) 



2A/a 

e 



^A 

e 



1 



BA 



F3*(0) 



2Mi 



Ga/3(0) 



2 

eA + i^2*(0) 



(^^3*(0) + F4*(0)) 



(8a) 

(8b) 

(8c) 

(8d) 



In the following, we will also use the relations that express the form factors 
^1 2 3 4 ill terms of the multipole form factors : 



R* = 



1 + r 
2 



Geo — -^tGe2 



Geo — -^tGe2 — 



Gmi — -rTGM3 

5 
4 

Gmi — -ttGmz 
5 



(1 



f: = - 



(1 



Geo — 
Geo - 



) Ge2 



gT I Ge2 



Gmi — ( 1 
Gmi — ( 1 



4^ 



G 



M3 



(9) 



For completeness we also express the form factors in terms of the covariant 
vertex functions ai,a2,ci and C2, used for instance in references [l2, 14. 11| : 



-C2 



-Fl = fll + a2 , F2 = -02 , = Cl + C2 , F4 

These relations, together with the identity 

u^ip'^Xya^'q^upip, A) = u^{p',X') (2Ma7^ - [p + pT)Mp^^) 



(10) 



(11) 



bring Eq. ([T]) into the form introduced in [ij, llJ, lU | • 

The empirical knowledge of the A electromagnetic moments is scarce, even 
though there were several attempts to measure the magnetic moment. The 
current Particle Data Group value of the magnetic dipole moment is given 
as [H: 

MA+ = 2.7+J3(stat.) ± 1.5(syst.) ± 3(thcor.)/xjv , (12) 

where fiN' = e/2A/jv is the nuclear magneton. This result was obtained from 
radiative photoproduction ('-/N — > ttN^') of neutral pions in the A(1232) region 
by the TAPS Collaboration at MAMI [16|, using a phenomenological model of 
the 7p — *■ tt'^p^' reaction [i3|. For the A+, Eq. (J2]) imphes : 

Gmi(0) = 3.5+J-^(stat.) ± 2.0(syst.) ± 3.9(theor.) . 



(13) 



The size of the error-bar is rather large due to both experimental and theoretical 
uncertainties. Recently a dedicated experimental effort is underway at MAMI 
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using the Crystal Ball detector [18[ , aiming at improving on the statistics of the 
TAPS data by almost two orders of magnitude. 

For the A electric quadrupole moment, no direct measurements exist. How- 
ever the electric quadrupole moment for the iV — > A transition has been mea- 
sured accurately from the ttN reaction at the A resonance energy [l^ : 

Qp^A+ = - (0.0846 ± 0.0033) c-fm^. (14) 

In the large- A'c limit of QCD, the two are related as follows [20| : 

'-^^0(1,]. (15, 



which, using the empirical value for Qp^^+ yields for the A+ quadrupole mo- 
ment: 



(0.048 ±0.002) e•fm^ (16) 



accurate up to corrections of order Note that using Eq. (j8dp . the value 

in Eq. HH) implies Ge2{0) = -1.87 ± 0.08. 

3. The 'natural' values of the e.m. moments 

As first argued by Weinberg [2l| , based on the Gerasimov-DrcU-Hearn sum 
rule, there is a 'natural' value for the magnetic moment of a pointlikc particle 
with spin, which corresponds to a gyromagnetic ratio g equal to 2. It has later 
been observed that all consistent field theories of charged particles with spin 



respect this value, see e.g. j22l.l23j. Given this universal result for the magnetic 
moment, it is reasonable to expect that all electromagnetic moments of pointlike 
particles are fixed at 'natural' values. To determine the values for the spin-3/2 
case, we examine the values of the e.m. moments of the gravitino in extended 
supergravity 24. 25l|. 



The gravitino, if it existed, would be a spin-3/2 particle described by a 
Rarita-Schwinger field which, in the framework of A/" = 2 supergravity, couples 
consistently to electromagnetism. Wc therefore expect that all the e.m. moments 
arising in this theory are 'natural'. To find their values, we start from the 
following Lagrangian density^ 

C = ViM7'^''"(i9o-e^„)V^,-m7^^7^>, 

+ em~^^p^{iniF^''' ~ K2J5F^"' )^P,. (17) 

It describes the spin-3/2 Rarita-Schwingcr field (V'/j) with mass m coupled to 
the electromagnetic field {A^j^) via the minimal coupling (with positive charge 



^In our conventions: 7'^'" = |[7'',7''], 7^"^" = §{7'"', 7"}, Fi^" = 9^ A'' - d^At^, Fi"" 
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e) and two non- minimal couplings ki and K2- As shown in Ref. [2y|, this is 
the most general Lagrangian of this type that gives the right number of spin 
degrees of freedom for a spin-3/2 particle. This theory, however, would still 



lead to rather subtle pathologies such as non-causal wave propagation [27|, |28[ , 
at least if no other fields are present. Adding gravity in a supersymmetric way 
makes the theory fully consistent from this viewpoint, but also constrains the 
non- minimal couplings as follows: 

Kl=K2 = 1. (18) 

We shall refer to these values as the 'SUGRA choice'. 
The e.m. vertex stemming from Eq. (|17|) is 



™ m 



where q = p'—p, and £0123 = +1- It is easy to verify that this coupling conserves 
the e.m. current: 

q^ua{p')r''''^{p\p)u^{p) = 0, (20) 
as well as, for the SUGRA choice, the supersymmetric current: 

(p'„-im7,)r"^^(p',p)u0(p)e^(g) = 0, (21a) 
Ua(p')r"''''(p',P)(P/3-|"^7/3)eM(9) = 0, (21b) 

where £^((7) is the photon polarization vector, and q ■ e = = in Eqs. (|2ip . 

The matrix elements of this vertex can be compared with the general de- 
composition of the spin-3/2 e.m. current given in Eq. ([T|), and in doing so we 
obtain the following result: 

F* = l + 2(Ki-HK2)T'"i'''l + 4r, (22a) 

F; = 2^*1 'I'"" 2, (22b) 

^^3* = 4(ki + At2) 8 , (22c) 

F* = 0, (22d) 

Thus, the values of gravitino's e.m. moments in A/" = 2 supergravity are: 

Geo{0)^1, Gmi(0)-3, G'£2(0) = -3, Gm3(0) = -1. (23) 

We take these values as the 'natural' values characterizing a structureless spin- 
3/2 particle, and interpret any deviation from them as a signature of internal 
structure. 

4. The 7* A A light-front helicity amplitudes 

In the following, we consider the electromagnetic A ^ A transition when 
viewed from a light front moving towards the A. Equivalently, this corresponds 
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to a frame where the baryons have a large momentum-component along the 
z-axis chosen along the direction of P = (p -\- p')/2, where p {p') are the initial 
(final) baryon four-momenta. We indicate the baryon light-front -I- component 
by (defining a* = a" ± a'^). We can furthermore choose a symmetric frame 
where the virtual photon four-momentum q has g"*" = 0, and has a transverse 
component (lying in the ly-plane) indicated by the transverse vector q±, sat- 
isfying q'^ = —q]_ = —Q^- In such a symmetric frame, the virtual photon only 
couples to forward moving partons and the -I- component of the electromagnetic 
current has the interpretation of the quark charge density operator. It is 
given by : J+(0) = -h2/3 m(0)7+u(0) - 1/3 (I(0)7+d(0), considering only u and d 
quarks. Each term in the expression is a positive operator since qj'^q oc |7'''(7p. 
We start by expressing the matrix elements of the ^"""(O) operator in the A 

as : 



(P+,^,AV+(0)|P+,-^,A) = (2P+)e*(^-^')^'AvA(Q'), (24) 



where A, A' denotes the A light-front helicities. and where q± = Q{cos(j)qex + 
sm(j)qey). The hclicity form factors Ay \ depend on only and can equiva- 
lently be expressed in terms of P* 2 3 4 '■ 



A3 3 
2 2 



^3 1 

2 2 



A3_l 
2 2 



^3_3 
2 2 



All 

2 2 



Ai_i 

2 2 



A_3_3=F* --F*, 

t1/2 

= A_l_3 ~ — Al3 = — — 
22 22 22 ./O 



A_3 1 = 3 = Al_3 = — 

2 2 2 2 2 2 /O 



-2P 



-A_ 3 3 



r3/2 FT 



A 1 , = {1--t]f* + - 



2P;' - P^ - r P 



Fl 



2^^F*+tF: 



- ( - - 2r ) P3* - 2rP« 



-1/2 r 



-A_ii 

2 2 



4Pi* - 2 (1 - 2t) P; - 2TP3* + T 



2r P; 



(25) 



Since the A ^ A electromagnetic transition is described by four independent 
form factors, one finds two angular conditions among the hclicity form factors 
of Eq. ^ : 

= (1 + At) ViA3 3- d>T^I'^ A3 i+2A3_i- ViAii, 
22 22 22 22 

1 3 

= 4x^/2^3 3 + \/3(l - 2T)yl3i H- -^3_3 - (26) 

22 22 222 922 
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5. The transverse charge densities for a spin-3/2 particle 

Wc define a quark charge density for a spin-3/2 particle, such as the A(1232), 
in a state of definite light-cone helicity A, by the Fourier transform : 

= r ^QMQb)A^^{Q^). (27) 

The two independent quark charge densities for a spin-3/2 state of definite 
helicity are given by pf{b) and pfib). Note that for a pointlike particle, the 

2 2 

'natural' values of Eq. lead to As 3 (Q^) = 1, implying 

Plib) = S^{b). (28) 

The above charge densities provide us with two combinations of the four inde- 
pendent A FFs. To get information from the other FFs, we consider the charge 
densities in a spin-3/2 state with transverse spin. We denote this transverse 
polarization direction by S± = coscjis^x + sin^sSy, and the A spin projection 
along the direction of S± by s±. We first express the transverse spin basis in 
terms of the helicity basis for spin-3/2. For the states of transverse spin s± = | 
and s_L = ^ this yields : 

-e'*-|A = -i)-^/3e2^^-|A = -^)|, (29) 

where the states on the rhs are the spin-3/2 helicity eigenstates. 

We can then define the charge densities in a spin-3/2 state with transverse 
spin s± as : 

p^^Jb) ^ |0^e-^^-^^^(P+,^,.^|J+(O)|P+,-^,.^).(3O) 
By working out the Fourier transform in Eq. (j30p for the two cases where s± = | 
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and 



i, using the A helicity form factors of Eq. (pS)) . one obtains 



and 



2ti 



+ 00 



dQ 
27r 



Q 



MQb)- M33 +3A11 
4V 22 22 

sin(0fc - (l)s) JiiQb) i (2V3A31 + 3Ai_i 

COs[2(0fc-05)] J2(Q6)^A3_i 

sin[3(0b-05)] J3(Q6) ^^3.3 



JoiQb)-[3As.+A..) 
sin(06-05) J (2\/3A3i 



(31) 



A, 



sin[3(0b-0s)] J3(Q&)^^|-| 



(32) 



where we defined the angle (pb in the transverse plane as, b = b{coB4)bex + 
sin^ijCj,). One notices from Eqs. (|31l32p that the transverse charge densities 
display monopolc, dipole, quadrupolc, and octupole field patterns, which are 
determined by the helicity form factors with zero, one, two, or three units of 
helicity flip respectively between the initial and final A states. 

It is instructive to evaluate the electric dipole moment (EDM) corresponding 
to the transverse charge densities ^ , which is defined as : 



7A - 



d%bp^,(b). 



Eqs. ([3TT321) yield : 



df = 3df = -(S^xeA {Gmi{0) - 3eA} 



\2Ma 



(33) 



(34) 



Expressing the spin-3/2 magnetic moment in terms of the (7-factor gives Gmi (0) = 
S'feA, so that the induced EDM df^ is proportional to 5 — 2. The same result 
was found before for the case of spin-1/2 particles in and spin-1 particles in 
Ref. 01 ■ One thus observes as a universal feature that for a particle without 
internal structure (corresponding with g = 2 [12, ll^), there is no induced EDM. 

We next evaluate the electric quadrupole moment corresponding to the trans- 
verse charge densities ■ Choosing S± = Bx, the electric quadrupole moment 
can be defined as : 



Qi 



d%{bl 



(&). 



(35) 
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From Eqs. (|31l32p one obtains 



Q|=-Qf = 2 {2[GMi(0)-3eA] + [G£;2(0) + 3eA]} (^^j . (36) 

We note that for a spin-3/2 particle without internal structure, for which Gi\ii{Q) = 
3eA and G£;2(0) = — 3eA according to Eq. ((23|) . the quadrupole moment of 
the transverse charge densities vanishes. It is thus interesting to observe from 
Eq. (|36p that is only sensitive to the anomalous parts of the spin-3/2 

magnetic dipole and electric quadrupole moments, and vanishes for a particle 
without internal structure. The same observation was made for the case of a 
spin-1 particle in Ref. 

At this point, it is important to emphasize the difference between densities 
and quadrupole moments defined in the lab and infinite momentum frames, and 
to discuss qualitatively the difference between Eq. ([5c)) and Eq. ([55)1 . For QCD 
in the lab frame, vacuum quark and gluon condensates as well as backward 
quark propagators play important roles in hadron structure. In contrast, the 
infinite momentum frame has neither condensates nor backward going propaga- 
tors Feynman diagrams. To boost a state from the lab to the infinite momentum 
frame is as complicated as solving for the state in cither frame, so at best, the 
effect of boosts can only be described in models. Because deep inelastic scatter- 
ing naturally specifies quark distributions and transverse densities in the infinite 
momentum frame, and because of the wealth of experimental information about 
parton distributions and generalized parton distributions, it is clearly valuable 
to develop a physical understanding of and intuition for distributions and den- 
sities in this frame, and the extraction of transverse densities from form factors 
can contribute substantially to this understanding and intuition. However, at 
the same time, it is essential to treat densities in this frame as physically distinct 
from those in the lab. 

In this context, it is useful to discuss three differences between the conven- 
tional quadrupole moment Qa of the three-dimensional density in the lab frame, 
Eq. (jScp . and the quadrupole moment Qf^ in the infinite momentum frame, 
Eq. l[36|) . The factor [Ga/i(0) — Sga] is an electric quadrupole moment induced 
in the moving frame due to the magnetic dipole moment, which is obviously 
absent in the lab. As discussed previously, the combination [Ge2{0) + 3e^] in 
Eq. p6p arises in the infinite momentum frame to describe non-point-like struc- 
ture because the value for a point particle is — 3eA, whereas no such term is 
pertinent in the lab. 

Finally, the overall factor 1/2 in Eq. ([36]) can be understood qualitatively by 
ignoring the distinction between the densities in the lab and infinite momentum 
frame, and considering a hypothetical limit in which one has a 3-dimensional 
density that is axially symmetric around the x axis and the transverse two- 
dimensional density is just the 2-dimensional projection of this density, without 
any change of physics between the lab and infinite momentum frames. Taking 
the spin axis along the x axis, the quadrupole moment Qsd of the 3-dimcnsional 
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distribution is defined as : 

Qzd = [ dxdydz{3x'^ - r'^) p3d{x,y,z), 

dxdydz [{x^ - y^) + {x^ - z^)] P3d{x, y, z). (37) 

For a 3-dimensional eharge distribution that is invariant under rotations around 
the axis of the spin, the two terms proportional to (x^ — y"^) and (x^ — z^) in 
Eq. ([37)1 give equal contributions yielding : 



Qid = 2 J dxdydz{x'^ ~ y'^) p3d{x,y,z). (38) 

Introducing the 2-diniensional charge density in the x'y-plane as : 

P2d{x,y) = J dzp3d{x,y,z), (39) 

one immediately obtains the relation 

Q3d^2Q2d, (40) 
with the quadrupole moment of the 2-dimensional charge density defined as : 

Q2d = dxdy {x^ ~ y^) p2d{x,y). (41) 



Because Q^d is proportional to GE2i^) in our hypothetical case, we see that 
Eq. (PP]) yields a Q2d that is half the value of Ge2{0), consistent with Eq. ([55]) . 
We will return to the differences between physics in the lab and infinite momen- 
tum frames and between Eq. ([8c]) and Eq. p6p in the discussion of models in 
the final section. 

In the same way as for the quadrupole moment, we can also evaluate the 
electric octupole moment corresponding with the transverse charge densities 
Pts±- Choosing S± = e^, the electric octupole moment can be defined as : 



Ot = ejd^b^sini3^,)p^,jb), 

= ej d'bby{3bl~bl)p^,Jb). (42) 



From Eqs. (j31l32p one obtains : 



2 32 



l{-GMi{0)-GE2{0) + GM3{0) + eA} (^^) (43) 



We note that for a spin-3/2 particle without internal structure, for which Gi\ii{Q) ■ 
3eA, Ge2{0) = — 3eA, and G7\/3(0) = — ca according to Eq. ((23)l . the electric 
octupole moment of the transverse charge densities vanishes. 
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6. Lattice calculation for the A (1232) e.m. form factors 



Since the A(1232) c.m. FFs are not known experimentally, apart from the 
scarce phenomenological information described in Section [21 we will rely on re- 
cent lattice QCD calculations [1^ for these FFs. For the A magnetic dipole mo- 
ment, first dynamical results, using a background field method, with Np = 2 + 1 
quark flavors were presented in Ref. [s^. For the A e.m. FFs, first dynamical 



results were presented in Ref. [ll|. We will discuss the latter calculations in 



more detail in the present section. 

We use three types of simulations that we refer to as SIM-I, SIM-II and 
SIM-III. In the SIM-I simulations, we use Wilson fcrmions and the standard 
Wilson plaquette gauge-action in the quenched approximation. In the SIM-II 
simulations, we use two dynamical degenerate flavors of Wilson fermions and 
the standard Wilson plaquette gauge-action. Finally, in the SIM-III simulations, 
we use the following mixed action scheme. For the dynamical sea quarks, we 
use the staggered Asqtad action by the MILC collaboration [3l[ and include two 
degenerate flavors of light quarks plus strange quarks. The strange quark mass is 
fixed to the physical strange quark mass. For the valence quarks we use domain 
wall fcrmions (DWF) that preserve a form of chiral symmetry on the lattice. 
The domain wall quark mass takes the values given in Tables [1] and has been 
tuned by adjusting the lightest pseudoscalar meson in the Asqtad calculation 31 1 
to have the same mass as the pseudoscalar meson using domain-wall fermions. 
Technical details of this tuning procedure are given in Refs. [s^, [s^. Using 
this mixed action, a number of interesting observables were recently evaluated, 
including the nucleon axial charge, where hybrid action chiral extrapolation to 
the physical point has been carried out [34| and moments of generalized form 
factors, which are also chirally extrapolated [l^- 

In this work we present a comparison between results with dynamical fermions 
obtained using the Wilson and the domain wall fermion discretization schemes. 
For finite lattice spacings these are different actions and agreement between 
them provides a non-trivial check of lattice cutoff effects. Results in the quenched 
approximation have the advantage that they are easily produced and have 
smaller statistical errors. Therefore they can be used to optimize the construc- 
tion of the three point functions. In addition, comparison with the unquenched 
data at a similar pion mass provides a measure of unquenching effects. For all 
simulations considered in this work the A is a stable particle. On the lattice we 
calculate the form factors Geq,Ge2, Gmi and Gm^ of which the electric charge, 
Gboj and the magnetic dipole, Gmi, form factors are dominant. 



6.1. Ground state dominance 

We use an interpolating field with the quantum numbers of the A^ baryon 
to create an initial state from the vacuum with non-zero overlap with the A+- 
state. The standard interpolating field that reduces to the quark-model wave 
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function in the non-rclativistic limit is given by 



^ abc 

V3 



(44) 

where C is the charge conjugation matrix. We note that this interpolating 
field creates both a spin-3/2 and a spin-1/2 state. However the overlap of this 
interpolating field with the spin-1/2 A+ that is higher in energy is very much 
suppressed and has negligible contribution to the correlator functions [35| . The 
overlap of the interpolating field with the spin-3/2 A+ is 

(f)|x<.a(0)|A(p,s)) =Zti,„(p,s), {Aip,s)\x.a{Om) ^Z*u^^ip,s), 

(45) 

and as in Eq. ([1]) denotes a Rarita-Schwinger spinor. Every vector component 
of this vector-spinor solves the free Dirac equation 

[^-Ma] w.(p,s) = 0, (46) 

and in addition the two auxiliary conditions 

p^u''{p,s) = and 7^u"(p,s) = (47) 

are satisfied. The well known expression for the spin-sum is given by 

3/2 

Ao-r = ^ Ua{p,s)Ur{p,s) 

s=-3/2 

_ f + Ma f Jcrlr 2p^Pt , Pair - Prla 



2MA V""'^ 3 SM'i UIa 



(48) 



using the normalization u°'Ua = — 1. 

To facilitate ground-state dominance, we employ a covariant Gaussian smear- 
ing [sl] on the quark-fields entering Eq. (|^^ 

q^it,x) = J2i^ + aH{x,y;U)r qp{t,y) (49) 
y 

3 

H{x,y;U) = ^ (C/^(a-,t)<5s,j;_^ + C/^(f ~ a^O-^.-^+a) (50) 

Here q is the local u or d quark field, q is the smeared quark field and [/^ is the 
SU (3)-gauge field. For the lattice spacings and pion masses considered in this 
work, the values a = 4.0 and n = 50 ensure ground state dominance with the 
shortest time evolution that could be achieved. 



6. 2. Correlation functions 

While the Breit frame is suitable for the introduction of multipolc form 
factors, and the infinite momentum frame allows for a clean definition of trans- 
verse charge densities, wc find that the most convenient kincmatical setup for 
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the lattice calculation is one where the final A-state is at rest {pf = 0). Fur- 
thermore, since lattice calculations are carried out in a Euclidean space-time, 
for the remainder of this chapter, all expressions are given with Euclidean con- 
ventions 



37|. 



We consider two-point and three-point functions constructed from the inter- 
polating fields defined in Eq. (|^ , 



= ^e--/po„ (xao(t,f/)w(o,o)) , (5i) 

G,^.(^^9,^) = ^e--?~r:;,„ (x.o(i/,x/)F^(i,x)x™'(0,0)) . (52) 



With lattice Wilson fcrmions, is the symmetrized, conserved electromagnetic 
current given by 

M^) = 1 ir.i^) + - '^)) - ^ Om(^) + - A)) where (53) 

3l[x) = g(a; + /i)i[7^ + l][/-i(a;)<7(a:)-fq(x)i[7^-l]t/^(a;)(z(x-K/i).(54) 

In the hybrid calculation we use the local vector current and determine the 
renormalization constant Zy by the condition Ge(^^ = 1 (in units of the electric 
charge) for the A+. The F matrices are given by 

r^ = J(l+7'), F'= = ^(l+7^)757fc, fc = l,2,3. (55) 

By inserting into the correlation functions complete sets of energy momentum 
eigenstatcs 

with ^ denoting all other quantum numbers like spin and its projection on some 
axis, one finds that the leading contributions for large Euclidean times t and 
tf — t are 

G^^(F^p,^) = ^^|Z|2e-^^('')*tr[F'^Af^(p)l + excited states, (57) 

G,^.(F^g,^) = |Zpe-*^-(*^-*)e-^-<-)*tr[F''Af,,(P/)Ofv.'Af..(p.)] 
-I- excited states . (58) 

To calculate the spin traces, the Euclidean versions of Eq. ^ and Eq. (|48p are 
employed. The unknown overlap Z-factors and the leading time dependence 
cancel in certain ratios of three- and two- point functions. Our preferred choice 
is the ratio 



Gfcfe(r4,0,</) V Gfefc(F4,0,</-i)Gfcfc(F4,pi,t)Gfcfc(F4,K,t/) ' 

(59) 
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with implicit summations over the indices k with k = 1,2,3. It becomes time 
independent for large Euclidean time separations tf — t and t: 



with 



2E 



A (p.) 



Ma 



Ml 



Ml 



Mi 



(60) 



(61) 



stemming from the state normalization and the 2-point function traces. 

Since we are evalua ting the correlator of Eq. (j52p using sequential inver- 
sions through the sink [38| . a separate set of inversions is necessary for every 
choice of vector and Dirac-indices. The total of 256 combinations is beyond our 
computational resources, and hence wc concentrate on a few carefully chosen 
combinations given below: 

n^H?) - E ^.^^n,/^,(r4,g) 

5iiEA + MA)C 



Gmi 



18Ml 



[Si.p.{q3 - 92) + ^2,m(9i - 93) + (53,^(92 - qi)] , 

(62) 



{Ea + 2Ma)C 



k=l 

—Geo 



3Ml 



[{Ma + EA)5i,t. + i9p(l - h.j.)] 



Ge2 ^^^o/i^^'^ [(Ma + Ea)^,^. + iq^{\ - <S4,^)] , 



(63) 



^f{q) = 



= G 



E ^jki^iiy^kij^^ ,q) 
j,k,i=i 

-iC 



E2 



3MI{Ea + Ma) 



Gmi 



Gm3 



G 



{qiq2 + 9293 + qsqi) 

X [(Ma + ^a)'54,m + 49^(1 - 
3 



6A/i(i?A + MA)f^^ 



E '^'^./^ 919293 ( 2 - 



91 + 92 + 93 - qk 
qk 



G 



30MI{Ea+Ma) 



k=l 



{16Ea + 14Ma)9i9293 



-10Ma(9i92 + 9293 + 939i)9fc 
'{8Ea + 7A/a)^^(9i + 92 + 93 - 9fe) 



9fe 



(64) 
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As expected, current conservation q^Hfj, = is manifest in the right hand side 
of the equations. From these expressions, all the multipole form factors can be 
extracted. For instance, Eq. (p^ is proportional to Gmi, while Eq. isolates 
Ge2 for /i = 4. Furthermore, these combinations are optimal in the sense 
that all momentum directions, each of which is statistically different, contribute 
to a given Q^-value. This symmetric construction yields a better estimator 
for the A-matrix elements than methods where only one momentum-vector is 
accessible. 



Wilson fermions 


V 


# confs 


(GeV) (GeV) 


(GeV) 




SIM-I 


: Quenched, P = 6.0, a"' = 2.14(6) GeV 




32'' X 64 
32^ X 64 
32^ X 64 


200 
200 
200 


0.1554 0.563(4) 0.645(9) 1.267(11) 
0.1558 0.490(4) 0.587(12) 1.190(13) 
0.1562 0.411(4) 0.503(23) 1.109(13) 


1.470(15) 
1.425(16) 
1.382(19) 


SIM-II: Unquenched, f3 = 5.6, a^' = 2.56(10) GeV 


24^ X 40 
24^ X 40 
24^ X 32 


185 [39| 
157 [39| 
200 [40| 


0.1575 0.691(8) 0.701(9) 1.485(18) 
0.1580 0.509(8) 0.566(12) 1.280(26) 
0.15825 0.384(8) 0.453(27) 1.083(18) 


1.687(15) 
1.559(19) 
1.395(18) 






SIM-III: Mixed action 




Asqtad 


(am„.d/3 = 


0.01/0.05), DWF (am„,<i = 0.0138), a'^ = 1.58(3) GeV 


V 


# confs 


(GeV) (GeV) 


niA 
(GeV) 


28^ X 64 


300 


0.353(2) 0.368(8) 1.191(19) 


1.533(27) 



Tabic 1: Parameters used in the calculation of the form factors. 



6.3. Sink-source separation 

Apart from the simulations themselves, sequential inversions are the most 
time consuming part of the calculation. They are done by constructing an 
appropriate sequential source, which requires that we fix the initial and final 
hadron states as well as the sink-source time separation, tf. Changing requires 
a new set of sequential inversions. Therefore we must first determine the optimal 
value of </ . The criterion is to choose i/ as small as possible so that statistical 
errors due to the exponential decrease of the signal are minimized but large 
enough so that we ensure that excited states with A quantum numbers are 
suppressed when the photon couples to a quark at time t from the source. 
Smearing techniques are essential for attaining early ground state dominance 
as discussed in the previous subsection. The optimal source-sink separation 
for a very similar calculation to ours was determined in Refs. [13, [l^l- The 
authors find that a distance of about 1 fm is sufficient to suppress excited state 
contributions well below our targeted statistical errors. This translates into a 
separation of 12 time shces for SIM-I, 13 for SIM-II and 8 for SIM-III. That this 
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Figure 2: The ratio defined in Eq. 1591 1 for a few of the lowest momenta. The left panel shows 
the combination from which Geo is extracted, whereas the right shows the one that isolates 
Gmi- This is for the quenched case where due to the high statistical accuracy excited state 
contributions are most visible. For higher momenta, as well as for the suppressed form factors 
the ratio of systematic to statistical error is lower. 



is indeed sufReient, is demonstrated in Fig. [2] and Fig. [3] where we show a few 
representative ratios with plateaus for SIM-I and SIM-III respectively. 

6.4- Volume dependence 

Another potential source of systematic error is the spatial size of our lattices. 
Given that for the quenched case we use a lattice of spatial size about 3 fm, 
we expect finite volume effects to be negligible. Similarly, for the mixed action, 
the spatial lattice size is = 3.5 fm giving Lgm^ = 6.4 ensuring small finite 
volume effects. In Ref. [l^] it was in fact shown that for the N to A matrix 
elements, finite volume effects are small for LsTUt^ ~ 4.5. Except for the lightest 
mass for dynamical Wilson fermions for which Lgm.,, ^ 3.6, for all our quark 
masses we have LgmT, > 4.6, and therefore we expect finite volume effects to be 
small. 

6.5. Data analysis 

For a given value of Q^, the combinations given in Eqs. ([5^ to (|64p are 
evaluated for all different directions of q resulting in this Q^, as well as for all 
four directions /.t of the current. This leads to an over-constrained linear system 
of equations, which is then solved in the least-squares sense yielding estimates 
of G^o, Ge2, Gmi and Gm3. This estimation is embedded into a jackknife 
binning procedure, thus providing statistical errors for the form factors that 
take all correlation and autocorrelation effects into account. 



18 



H=4, p=(0,0,0)- 
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Figure 3: The same as in Fig. [2] but for the hybrid case. 



6. 6. Simulation parameters 

The simulation parameters of our calculation are given in Table [1] In the 
case of Wilson fcrmions. the lattice spacing has been estimated using the nu- 



cleon mass in the chiral limit 41|. The value extracted using the nucleon mass 



is in agreement with the determination using the hadronic scale rg defined via 
the force between static quarks at intermediate distance [i^. The lattice spac- 
ing in the case of staggered fermions has been determined from heavy-quark 
spectroscopy |43| as a = 0.1241 fm with a statistical uncertainty of 2%. The 
calculation has been performed on relatively large lattices in the quenched ap- 
proximation and in the case of the mixed action of spatial extent 2.9 fm and 
3.4 fm respectively. On the other hand, in the case of dynamical Wilson fermions 
we use a finer lattice with a spatial size of only 1.9 fm. 

Isospin symmetry relates results obtained for the to those for the A^"*" , A° 
and A~ since they differ only by a charge-factor. The pion and A masses cor- 
responding to the three values of the hopping parameter, k, considered in sim- 
ulations SIM-I and SIM-II as well as the smallest available pion mass in the 
mixed action are summarized in Table [TJ The results presented in the next 
subsection are obtained from the connected contributions to the full electro- 
magnetic current. These yield the iso- vector form factors with the iso- vector 
current being given by = ^7^71 — d'^^j^d. In this case the disconnected dia- 
grams cancel. The iso-scalar form factors have contributions from disconnected 
diagrams. At present, it is not yet computationally feasible to calculate the con- 
tributions arising from disconnected diagrams since they require the calculation 
of all-to-all propagators. An exact evaluation of all-to-all propagators requires 
spatial volume inversions for each time slice t that we are interested in. Such 
a computation is not feasible. Stochastic techniques have been developed to 
estimate all-to-all propagators [4^ . Various improvements are being developed 
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Figure 4: Comparison of three different QCD lattice calculations for the A+(1232) form 
factor Geo- The upper curve is for SIM-II and SIM-III using dynamical quarks, whereas the 
lower curve is for SIM-I using the quenched approximation 12911 . The lines show the fits to a 
dipole form of the lattice results at the smallest pion mass. For the case of dynamical quarks 
(upper graph) the fit is made to the results obtained using the mixed action. The error band 
is calculated using a jackknife analysis of the fit parameters. 



in order to improve the measurements [45j . However, these disconnected con- 
tributions are particularly hard to calculate not just because they require the 
all-to-all propagators but also because they are noise-dominated [i^. A cal- 
culation of disconnected contributions in the case of the electromagnetic form 
factors have shown that these are consistent with zero within large statistical 
errors 47|. We expect that, like in the case of nucleon electromagnetic form 
factors, the disconnected contributions to the electromagnetic A form factors 
are small. Therefore assuming that the disconnected diagrams are small then 
we obtain the electromgnetic form factors of A. 

6.7. Results 

We show in Figs. |4] and [5] the two dominant Sachs form factors Geo and 
Gmi as a function of the for all the light quark masses. In Fig. [6] we display 
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Figure 5: Comparison of three different QCD lattice calculations for the A+(1232) magnetic 
dipole form factors Ga/i- The lines show the fits to an exponential form of the lattice results 
at the smallest pion mass. The rest of the notation is the same as that in Fig. |3] 



the Q^-dependence of the subdominant Sachs form factor Ge2- Their values 
arc tabulated in Tables [3l [4] and [5] of the Appendix. 

For Geo, we consider a dipole parameterization of the lattice results : 



' - (TTwXi^' '''' 



with A|;q treated as a free parameter. The choice of the parameterization for the 
A multipole form factors given in Eq. (|65|) ensures that the helicity conserving 
form factors As 3 and Aii behave as l/Q^ for large Q^. For Gmi and Ge2, 
we instead consider exponential parameterizations with two free parameters for 
each since the expected large Q^-dependence for these form factors is stronger 
than for a dipole : 

Gmi{Q^) = GA/i(0)e-«'/^'", 

Ge2{Q^) = GE2{0)e-Q'/^'^^-. (66) 
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Figure 6: Comparison of three different QCD lattice calculations for the A+(1232) electric 
quadrupole form factors Ge2- The notation is the same as that in Fig. [5] 



The mag netic octupolc form factor Gms evaluated in the quenched approxima- 
tion [29[ is found to have a smaU negative value albeit with large error. The 
results for Gm3 using dynamical Wilson fermions obtained in this work are 
compatible with Gm3{Q^) — within the current statistical accuracy. We will 
therefore not discuss further lattice results for Gms- 

The values of the parameters entering the dipole fit of Eq. ((65)) for Geo, and 
the exponential fits of Eq. for Gmi and Ge2, fitted to the three different 
lattice calculations SIM-I, SIM-II and SIM-III are listed in Tabled 

In Figs. m Eland [5] we show the fits to the lattice results for the form factors 
Geo, Gmi and Ge2 for the smallest pion mass in each of the three type of 
simulations. One sees that for Geo all three calculations give similar results. 
For Gmi and more so for Ge2, ^ larger spread in the resulting fits is observed. 
For Ge2{0) the hybrid lattice calculation yields a quadrupole moment consistent 
with the large- A^c value, Eq. pB]) . whereas the quenched and dynamical Nf = 2 
Wilson calculations yield only about half this value. It should, however, be noted 
that the present dynamical calculations still have substantially larger error bars 
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rn^ 
(GeV) 


(GeV2) 


Ga/i(0) 


Ami 
(GeV^) 


Gb2(0) 


(GeV^) 


SIM-I 


0.563(4) 
0.490(4) 
0.411(4) 


1.243(27) 
1.173(28) 
1.109(29) 


2.695(59) 
2.678(72) 
2.635(94) 


1.015(32) 
0.994(36) 
0.975(23) 


-0.46(20) 
-0.52(26) 
-0.81(29) 


1.06(42) 
0.91(39) 
0.696(200) 


SIM-II 


0.691(8) 
0.509(8) 
0.384(8) 


1.667(106) 
1.318(48) 
1.144(54) 


2.589(78) 
2.68(13) 
2.35(16) 


1.300(81) 
0.977(74) 
1.015(76) 


-0.71(49) 
-1.68(88) 
-0.87(67) 


1.00(64) 
0.76(30) 
1.75(1.39) 


SIM-III 


0.353(2) 


1.160(78) 


3.04(24) 


0.935(122) 


Z.UD_2 35 


>S4+^ '^^ 



Table 2: Parameters for the A+(1232) FFs for three different lattice QCD calculations : 
quenched Wilson (SIM-I), dynamical Nj = 2 Wilson (SIM-II), and a hybrid calculation (SIM- 
III) described in the text. For Geo, the dipole parameterization of Eq. I|65|l is used. For Gmi 
and Gb2i the exponential parameterization of Eq. I|66p is used. The errors on the parameters 
are jackknife errors. 

than the quenched calculations. 

7. The A+ transverse densities 

In this section, we will use the values of the parameters determined in the 
lattice fits to the A e.m. FFs, discussed in Section HI to evaluate the quark 
transverse densities in the A+(1232). 

In Fig.[7l we compare the A"*" transverse densities in helicity states of A = 3/2 
and A = 1/2, using Eq. (|27p . for the quenched lattice calculations. A comparison 
reveals that both are very similar, showing a positive core in the transverse 
charge density of about 0.7 fm. One notices that the quark transverse density 
in a A = 1/2 state is slightly more concentrated than its A = 3/2 counterpart. 
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Figure 7: Quark transverse charge densities in a A+(1232) of definite Ught -cone helicity. 



Upper left panel 



Upper right panel : pj . The light (dark) regions correspond to the 



largest (smallest) values of the density. The lower panel compares the density along the j/-axis 
for (dashed curve) and p'^ (solid curve). For the A c.m. FFs, the quenched lattice QCD 

2 2 

results are used (fit of Table [2} . 



In Fig. [HI the transverse densities of Eqs. (j31|32p arc compared for a A+ 
that has a transverse spin. It is seen that the quark charge density in a A"*" in 
a state of transverse spin projection s± ~ +3/2 is elongated along the axis of 
the spin (prolate deformation) whereas in a state of transverse spin projection 
s_L = +1/2 it is elongated along the axis perpendicular to the spin. 
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Figure 8: Quark transverse charge densities in a A"*" (1232) which is polarized along the positive 
X-axis. Upper left panel : 3 . Upper right panel : j . The light (dark) regions correspond 

to the largest (smallest) values of the density. In order to see the deformation more clearly, a 
circle of radius 0.5 fm is drawn for comparison. The lower panel compares the density along 
the j/-axis for 3 (dashed curve) and p^ ^ (solid curve). For the A e.m. FFs, the quenched 

lattice QCD results are used (fit of Table [2]l. 



The corresponding dipole, quadrupolc and octupole field patterns in the 
transverse quark charge density for a transversely polarized A are shown in 
Fig.l 
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Figure 9: Multipole field patterns in the quark transverse charge density 3 in a A+(1232) 

that is polarized along the positive x-axis. Upper left panel : dipole field pattern; upper right 
panel : quadrupole field pattern; lower panel : octupole field pattern. For the A e.m. FFs, 
the quenched lattice QCD results are used (fit of Table |2j. 



To convey the consistency of all three types of simulations, we show the 
differences in the transverse densities between the quenched calculation and the 
two dynamical lattice calculations in Fig. [TOl One notices that the differences 
are very small and arise primarily in the central densities for quenched and 
dynamical Wilson calculations. 
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P3I2^ [1/fm'] 





Figure 10: Comparison of the densities along the j/-axis for (upper panel) and 3 (lower 

2 2 
panel), for three different QCD lattice calculations of the A+(1232) e.m. FFs according to 
the fit of Table [2] Solid (red) curves : quenched Wilson result; long dashed (blue) curves : 
dynamical Nf = 2 Wilson result; dotted (black) curves : hybrid lattice calculation. 

8. Discussion and Conclusions 

In this paper we have studied the electromagnetic properties of the A(1232)- 
resonance using lattice QCD. Lattice results for the A electromagnetic form 
factors have been presented for pion masses down to approximatively 350 MeV 
for three cases: quenched QCD, two flavors of dynamical Wilson quarks, and 
three flavors of quarks described by a mixed action combining domain wall 
valence quarks and staggered sea quarks. 
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We have also presented the formaUsm for understanding the spin-3/2 A form 
factors in terms of precisely defined transverse quark charge density distribu- 
tions in the infinite momentum, or light-front, frame. For reference, we showed 
that a point-like transverse charge density, that is, a two dimensional delta 
function, yields the specific values G_eo(0) = 1, Gmi(0) = 3, Ge2{0) = —3, 
and Gm3(0) = — 1 and that these values are the same as the 'natural' values 
occurring for structureless gravitinos in extended supcrgravity. This strongly 
supports the hypothesis that these are universal values characterizing struc- 
tureless spin-3/2 particles. The fact that our lattice calculations yield values 
significantly different from these 'natural' values indicates the presence of inter- 
esting, nontrivial structure, and we have used the helicity amplitudes extracted 
from the electromagnetic form factors to construct the two-dimensional trans- 
verse quark charge densities for longitudinally and transversely polarized spins. 
A salient feature is that for the state with transverse spin projection 3/2, the 
transverse quark charge distribution is prolately deformed along the axis of the 
spin. 

Since there is also a very different language for describing the internal de- 
formation of the A in terms of a three-dimensional intrinsic state density, it is 
useful to briefly summarize the status of model studies using that description 
and discuss the differences between the two languages. In nuclear physics, for 
which the nucleus is sufficiently massive that for momentum transfers charac- 
teristic of the spatial size of the nucleus, electromagnetic form factors are well 
described by the Fourier transform of a three-dimensional charge density, there 
is a fruitful history of determining the deformations of intrinsic states. The basic 
idea [i^ is to treat collective rotations of a well-deformed nucleus in terms of a 
deformed intrinsic state that is axially symmetric in a body fixed-frame in which 
the projection of the total angular momentum J along the body-fixed symme- 
try axis is K. Following Ref. [48| . the relation between the intrinsic quadrupole 
moment, Qo, and the spectroscopic quadrupole moment, Q, measured experi- 
mentally and corresponding to the Fourier transform of the density in the lab 
frame, is 

Q- (j + i)(2J + 3)^°- ^^^> 

In nuclear physics, this approximation is not only successful phcnomcnologically, 
but also ab initio nuclear mean-field theory calculations using many-body theory 
that accurately predict electromagnetic form factors of spherical nuclei also 
predict the form factors arising from deformed intrinsic states in rare-earth and 
actinide nuclei [iot . 

It is natural to ask whether this intrinsic state approximation is also applica- 



ble in studying the nucleon. One obvious example is the Skyrme model jSOl l5l| , 
in which a single intrinsic soliton state generates an infinite band of collective 
spin-isopin rotations, with the spin-1/2 isospin-1/2 nucleon being the lowest 
state and the spin-3/2 isospin-3/2 A being the first excited state. It is clear, 
however, that the Skryme model cannot be a quantitative model of the nucleon 
and A, since the infinite tower of predicted additional states are not observed 
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in Nature. Buchmann and Henley [52] reviewed a variety of models of the A, 
and concluded that for a quark model, a collective model, and a pion cloud 
model, the intrinsic state was oblate, corresponding to a negative value for Qg- 
Studying a sphcriodal bag model, VioUier, Chin, and Kerman (ssj concluded 
that the K = 3/2 intrinsic state (the same as considered by Buchmann and 
Henley) is oblate and that the K = 1/2 intrinsic state is prolate (as in the 
Skyrme Model). We note that for a J = 3/2 A, Eq. (|67)) specifies that for a 
if = 3/2 intrinsic state Q and Qo have the same sign whereas for a if = 1/2 
state they have opposite signs. Thus, both a prolate K = \/2 intrinsic state and 
oblate if = 3/2 state yield a negative spectroscopic quadrupole moment, and 
indeed, all the results in Refs. [s^, [H^ are consistent with a negative value for 
Q and, thus by Eq. ((8c)) . with a negative value for Ge2 as we have calculated in 
Table [2l A recent calculation of the A+ form factors in the chiral quark soliton 
model [SJ] found a value G_e2(0) = —2.145, consistent with the value of the 
hybrid lattice result shown in Table [21 

Unfortunately, it is not possible to go further and establish a correspondence 
between our result of a prolate deformation of the 2-dimensional transverse spin 
density along the axis of the spin-3/2 state and the model statement of an 
oblate intrinsic 3-dimensional density for a if = 3/2 intrinsic state. As we have 
emphasized, even a structureless A would have Ge2 = —3, so a negative Ge2 
less negative than —3 would indicate prolate deformation in our case while still 
being consistent with oblate deformation for an intrinsic state. It is possible 
that study of a simple relativistic model that could be boosted into the infinite 
momentum frame could go somewhat further in bridging the gap between these 
two different descriptions. 
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A. Lattice results for the Delta electromagnetic form factors 

In Tables mm and[5l we list the lattice results for the Delta e.m. FFs in the 
quenched; dynamical Wilson, and hybrid calculations respectively. 
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(GoV^) 


G EO 




G E2 




0.173093(69) 


0.7750(27) 


2.340(47) 






0.33976(26) 


0.6228(41) 


1.913(40) 


-0.37(13) 




0.50068(54) 


0.5136(54) 


1.600(36) 


-0.250(97) 




0.65640(89) 


0.4368(68) 


1.361(32) 






0.8074(13) 


0.3728(69) 


1.182(29) 


-0.216(72) 


m,r = 0.563(4) GeV 


0.9541(18) 


0.3205(75) 


1.027(29) 


-0.170(60) 




1.2358(28) 


0.2527(85) 


0.809(31) 


-0.177(57) 




1.3715(34) 


0.2174(87) 


0.701(31) 


-0.159(50) 




1.5040(40) 


0.194(10) 


0.628(33) 


-0.084(63) 




1.6335(46) 


0.1698(97) 


0.549(31) 


-0.080(50) 




1.7603(53) 


0.142(11) 


0.460(38) 


-0.136(57) 




1.8845(59) 


0.140(12) 


0.439(35) 


-0.100(55) 




2.0063(66) 


0.117(10) 


0.372(31) 


-0.085(44) 




2.2430(79) 


0.084(16) 


0.284(47) 






0.172859(82) 


0.7646(32) 


2.324(58) 






0.33890(30) 


0.6076(46) 


1.889(48) 


-0.41(16) 




0.49887(64) 


0.4965(61) 


1.575(43) 


-0.27(12) 




0.6534(11) 


0.4212(77) 


1.329(37) 






0.8030(15) 


0.3574(75) 


1.155(34) 


-0.209(85) 


= 0.490(4) GeV 


0.9481(21) 


0.3055(83) 


1.004(33) 


-0.160(70) 




1.2263(33) 


0.2428(90) 


0.797(35) 


-0.196(68) 




1.3600(40) 


0.2060(93) 


0.682(35) 


-0.167(61) 




1.4906(47) 


0.182(11) 


0.608(37) 


-0.054(76) 




1.6181(54) 


0.159(10) 


0.531(34) 


-0.062(59) 




1.7427(61) 


0.131(12) 


0.440(41) 


-0.151(72) 




1.8648(68) 


0.133(12) 


0.432(39) 


-0.118(65) 




1.9844(76) 


0.110(10) 


0.358(33) 


-0.091(50) 




2.2167(91) 


0.073(16) 


0.259(47) 






0.172616(98) 


0.7546(42) 


2.320(78) 






0.33800(36) 


0.5933(55) 


1.876(62) 


-0.50(22) 




0.49700(75) 


0.4797(72) 


1.558(54) 


-0.32(15) 




0.6503(12) 


0.4061(92) 


1.295(48) 






0.7985(18) 


0.3422(84) 


1.131(41) 


-0.20(12) 


= 0.311(4) GeV 


0.9420(24) 


0.2908(95) 


0.987(41) 


-0.170(87) 




1.2167(39) 


0.236(10) 


0.797(41) 


-0.244(90) 




1.3485(46) 


0.196(10) 


0.665(41) 


-0.199(83) 




1.4770(54) 


0.172(12) 


0.590(43) 


-0.03(10) 




1.6025(62) 


0.150(11) 


0.518(39) 


-0.041(79) 




1.7250(70) 


0.121(13) 


0.417(48) 


-0.19(11) 




1.8450(79) 


0.127(13) 


0.436(46) 


-0.193(88) 




1.9624(87) 


0.103(11) 


0.349(37) 


-0.117(62) 




2.190(10) 


0.066(17) 


0.231(52) 





Table 3: Delta form factors from quenched Wilson fcrmions. 
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(GoV^) 


G EO 


G]\J~l 


G E2 




0.43260(36) 


0.6248(40) 


1.900(50) 






0.8364(13) 


0.4352(65) 


1.348(44) 


-0.33(15) 




1.2164(25) 


0.337(11) 


0.998(49) 


-0.17(12) 




1.5765(40) 


0.282(19) 


0.761(56) 




m,r = 0.691(8) GeV 


1.9195(56) 


0.235(22) 


0.638(58) 


-0.12(10) 




2.2475(74) 


0.189(21) 


0.501(61) 


-0.095(79) 




2.866(11) 


0.125(32) 


0.312(84) 


-0.012(96) 




3.159(13) 


0.103(24) 


0.225(56) 


-0.043(75) 




3.442(15) 


0.069(19) 


0.136(43) 


-0.155(96) 




3.717(17) 


0.054(27) 


0.129(58) 


-0.04(11) 




3.984(19) 


0.081(47) 


0.106(70) 


-0.12(18) 




4.244(21) 


0.034(13) 


0.087(32) 


-0.008(58) 




4.497(23) 


0.040(33) 


0.078(72) 


-0.010(93) 




4.985(27) 


0.019(34) 


-0.004(49) 






0.42932(57) 


0.5748(85) 


1.743(76) 






0.8250(20) 


0.3785(82) 


1.199(55) 


-0.49(21) 




1.1939(39) 


0.288(13) 


0.745(60) 


-0.41(19) 




1.5409(61) 


0.241(30) 


0.546(78) 




m,r = 0.509(8) GeV 


1.8694(85) 


0.181(16) 


0.433(55) 


-0.15(12) 




2.182(11) 


0.194(51) 


0.36(11) 


-0.04(13) 




2.768(16) 


0.093(33) 


0.181(75) 


-0.22(14) 




3.044(19) 


0.037(18) 


0.043(37) 


0.004(91) 




3.311(22) 


0.019(15) 


0.039(26) 


-0.021(84) 




3.569(25) 


0.021(39) 


0.032(62) 


-0.06(14) 




3.819(27) 


0.045(41) 


-0.009(59) 


0.05(17) 




4.062(30) 


0.013(16) 


0.020(29) 


-0.018(74) 




4.298(33) 


0.003(19) 


-0.001(27) 


0.099(82) 




4.753(38) 


0.015(43) 


-0.00(11) 






0.42485(79) 


0.566(15) 


1.530(92) 






0.8099(26) 


0.368(16) 


1.048(62) 


-0.53(27) 




1.1647(50) 


0.276(32) 


0.891(97) 


-0.41(24) 


= 0.384(8) GeV 


1.4953(77) 


0.195(48) 


0.57(11) 






1.806(11) 


0.141(80) 


0.354(39) 


-0.30(15) 




2.100(14) 


0.11(12) 


0.337(61) 


-0.27(15) 




2.648(20) 


0.11(25) 


0.062(99) 


-0.42(49) 




2.905(23) 


0.05(31) 


0.095(37) 


0.04(12) 




3.152(26) 


0.04(39) 


0.084(53) 


-0.48(24) 




3.391(29) 


0.03(47) 


0.047(95) 


-0.19(39) 




3.622(32) 


0.04(57) 


-0.040(94) 


0.17(23) 




3.846(35) 


0.03(68) 


0.001(54) 


-0.05(20) 




4.063(38) 


0.04(78) 


0.014(35) 


-0.06(10) 




4.481(44) 


-0.0(1.0) 


0.011(69) 





Table 4: Delta form factors from dynamical Wilson fermions. 
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(GeV^) 


G EO 




G E2 


0.124094(51) 


0.8223(95) 


2.83(21) 




0.24512(19) 


0.686(13) 


2.26(17) 


-1.47(88) 


0.36330(42) 


0.584(17) 


1.99(17) 


-0.84(72) 


0.47881(71) 


0.499(22) 


1.89(17) 




0.5918(11) 


0.432(21) 


1.53(14) 


-0.76(48) 


0.7025(15) 


0.386(23) 


1.38(14) 


-0.56(46) 


0.9175(24) 


0.310(29) 


1.15(14) 


-0.20(44) 


1.0220(29) 


0.287(35) 


1.04(15) 


-0.47(43) 


rrin = 0.353(2) GeV 1.1246(35) 


0.239(31) 


0.83(13) 


-0.05(47) 


1.2255(41) 


0.223(34) 


0.78(14) 


-0.44(47) 


1.3247(47) 


0.261(88) 


0.90(31) 


-0.72(65) 


1.4223(53) 


0.187(35) 


0.67(12) 


-0.38(43) 


1.5184(59) 


0.176(45) 


0.65(16) 


-0.81(48) 


1.7064(73) 


0.094(42) 


0.27(15) 




1.7983(80) 


0.125(45) 


0.41(13) 


-0.64(48) 


1.8890(87) 


0.095(34) 


0.327(85) 


-0.40(39) 


1.9785(94) 


0.101(40) 


0.30(13) 


-0.84(64) 


2.067(10) 


0.073(35) 


0.276(87) 


-0.39(44) 


2.154(11) 


0.073(37) 


0.191(84) 


-0.51(43) 


2.240(12) 


0.124(89) 


0.23(27) 


-1.3(1.4) 


2.409(13) 


0.079(55) 


0.16(16) 


-0.45(71) 



Table 5: Delta form factors from the hybrid approach. 
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